Spectrum and thermal fluctuations of a microcavity polariton Bose-Einstein 
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The Hartree-Fock-Popov theory of interacting Bose particles is developed, for modeling exciton- 
polaritons in semiconductor microcavities undergoing Bose-Einstein condensation. A self-consistent 
treatment of the linear exciton-photon coupling and of the exciton non-linearity provides a thermal 
equilibrium description of the collective excitation spectrum, of the polariton energy shifts and of 
the phase diagram. Quantitative predictions support recent experimental findings. 
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A major advance in the research on quantum fluids 
was made with the recent experimental observation of 
quantum degeneracy and off-diagonal long-range order 
(ODLRO), in a gas of exciton-polaritons in a semicon- 
ductor microcavity [H, 0, [1] . Bose-Einstein condensation 
(BEC) is the most appealing way of interpreting these 
findings. However, the 2-D nature of the system, the 
presence of disorder, the hybrid exciton-photon composi- 
tion of polaritons, and the composite nature of excitons, 
call for models that account for the peculiar aspects of the 
polariton gas. Along this line, recent theoretical works 
have been successful in describing many specific aspects 
of the system. Disorder at the quantum well interfaces, 
in particular, was accounted for by a BCS-likc theory in 
which bound excitons were modeled as particle- hole ex- 
citations coupled to the photon field @, @, @] . For strong 
exciton disorder, this theory models the density depen- 
dence of the polariton spectrum and shows how the linear 
optical response can probe successfully the phase tran- 
sition. More generally, in presence of localization, the 
lower-energy region of the density of states is modified 
with respect to an ideal 2-D system 0, H| , making BEC of 
a trapped polariton gas the most suited description of the 
system. Finally, the short radiative lifetime of polaritons 
results in non-equilibrium effects In & [IS EI IS Ei ELI, 
leading to condensate depletion [14| and to a modified ex- 
citation spectrum 0, [l3| . These effects should however 
be weak at densities above the condensation threshold 
, and are expected to be negligible for a polariton life- 
time longer than fO ps.[24| Very recent measurements 
clearly show that thermal equilibrium polariton BEC can 
indeed be achieved 0] . In spite of the high relevance of 
the existing theoretical frameworks, a basic question re- 
mains still unanswered. Are the experimental findings 
correctly interpreted in terms of a quantum field the- 
ory of interacting bosons? Such a theory should account 
sclf-consistcntly for the linear coupling between two Bose 
fields - photons and excitons - and for the Coulomb and 
Pauli non-linearities arising from the composite nature of 
excitons. 

In this Letter, we answer this question by general- 
izing the Hartree-Fock-Popov (HFP) [H, [3| theory of 



BEC to the case of two coupled Bose fields. In order to 
address the fundamental thermodynamical properties of 
the polariton gas, we assume thermal equilibrium - as 
was done in previous works based on non-bosonic mod- 
els 0, [H[ - and discuss kinetic effects elsewhere Q • The 
system, including the Coulomb and Pauli non-linear exci- 
ton terms, is described within an effective boson Hamil- 
tonian [o, 0, valid well below the exciton Mott 
density. We derive coupled equations for the condensate 
wave function and the field of excitations, and study the 
solutions for parameters modeling a recent experiment 
0. We discuss the collective excitation spectrum, the 
density-dependent energy shifts, the onset of off-diagonal 
long-range order and the phase diagram. Our analysis 
provides a generally good account of the experimental 
findings, in particular by reproducing the measured crit- 
ical density and the energy shifts of the two polariton 
modes. 

Let us introduce the exciton and photon operators b k 
and Cfe, obeying Bose commutation rules. 25j The polari- 
ton Hamiltonian then reads 



H = if o + Hji + H x + H s 



(1) 



where H = ^ k £^j c t k + e^c^Ck is the non- 
interacting term, Hr = TiQr :X)k(^k^k + h.c.) de- 
scribes the linear exciton-photon coupling, and H x = 
lEfcfc' g u z( k > k ''q)^k+q^k'-q^k'&k is an effective 2-body 
exciton Hamiltonian, modeling both Coulomb interaction 
and the effect of Pauli exclusion on electrons and holes. 



The term H s = J2 kk > q w «( k > k '> q)(c k+q ^L- q 6 k'&k + h.c.) 
models the effect of Pauli exclusion on the exciton os- 
cillator strength fl7l , which decreases for increasing 
exciton density [211 ]. We account for the full momentum 
dependence of v x (k., k', q) and v s (k, k', q) 17, 3|. They 
vanish at large momenta, preventing the ultraviolet di- 
vergence typical of a contact potential 22|.[26| For clar- 
ity, however, we use the short form u x (k, k', q) — ► v x and 
v s (k, k',q) — > MlR/n s in what follows, where n s is the 



saturation density of the exciton oscillator strength [17 1 . 

As non-interacting exciton and photon modes, we as- 
sume states in a square box of area A with periodic 
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boundary conditions. In real microcavities, interface fluc- 
tuations affect the photon modes, resulting in a disorder 
potential and in polariton localization over a few pm • 
For equilibrium BEC, the energy spacings between the 
lowest-lying states determine locally the effect of quan- 
tum fluctuations. Our finite size assumption thus models 
in a simple way the local structure of the spectrum of the 
disordered system. 

Via the Bogolubov ansatz, the exciton and the photon 
fields are written as 



# x ,c(r, t) = $ K , c (r, t) + 7/>z,c(r, t), 



(2) 



i.e. as the sum of a classical symmetry-breaking term 
$> x>c (r, t) for the condensate wave function, and of a 
quantum fluctuation field ip XjC (r,t). Heisenberg equa- 
tions of motion result in two coupled equations for 
$ XiC (r, t). In the HFP approach, anomalous correlations 
of the excitation field are neglected HI, 1(| and we obtain 



2 Re {n x 



m R ( i _ ^ \ $ 

n. 



e£$ c + nn R [ 1 



J + v (r 



4> 



(3) 



Here, we introduce the density matrix n 

n X 6 (x,£ = x,c), 



n 



where 



<f>*$5 and 



Ek^o n x?( k ) = Ek#o(°x( k )°?( k )) arc thc contribu- 
tions from the condensate and from the excited states, re- 
spectively. In our notation, O^k) = Of(k) and 0|(k) = 

0|(— k), with O x = b, and O c = c. The quantities n x {(k) 
are computed self-consistently as described below. By 
setting $x,c(t) = e~ lEt/h $ X:C {Q) into (J3J, we obtain a 
generalized set of two coupled Gross-Pitaevskii equations 
for the condensate eigenstate. The two solutions of ([3]) 
correspond to the lower and upper polariton respectively, 



and can be expressed as $ 
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thus fixing thc Hopficld coefficients of thc polariton con- 
densate. The low-energy solution E l Q p = fi defines the 
chemical potential. 

In the HFP theory, Eqs. ([3]) are coupled to the 
field-equations for excitations. In analogy with the 
case of a single Bose gas, we define the 4x4 matrix 
G(k, iu) n ) = {g*i (k.,iuj n )}*f~i2 i whose elements are the 
thermal propagators of the excited particles [l6[ : 



3ji 



x? (k,^„) 



dre^(6U-k,T)dUk,0y) r , (4) 



where huj n = 2nn//3, n = 0, ±1, ... are the Matsubara en- 
ergies for bosons and (...} T represents the thermal average 
of the imaginary-time ordered product. 

The propagator matrix G(k, iu> n ) obeys the Dyson- 
Belaev equation 

G (k, iuj n ) = G° (k, iio n ) [1 + E (k, iuj n ) G (k, iu n )] , (5) 



where G° = {$ (k, iw n )}f = SM(-)Hu; n ~e^ + p}- 1 
is the matrix of the unperturbed propagators, and 



S(k, iu n ) 



Z* x (k,iu n ) Z xc (k,iu> n ) 
E CT (k,^„) E cc (k,iw„) 



(6) 



is the 4x4 self-energy matrix, here written in a 2 x 2-block 
form. In thc HFP limit, the self-energy is independent of 
frequency and reads 
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while T/jf = Yijf and E^ = 0. For each value of k, the 
analytic continuation of the Green's functions (k, z) 



have four poles at z = ±E l ^ up \ They are the posi- 
tive and negative eigen-energies of the lower- and upper- 
polariton Bogolubov modes. The components of the cor- 
responding eigenvectors h^(fe) = (X u ,C u ,X v ,C v ) l (J (k) 
(i = Ip, up and a = +, — ) are the generalized Hopficld co- 
efficients of thc normal (X u , C u ) and anomalous (X v , C v ) 
kind, with normalization iX^^-lX^p + ICup-IGt,! 2 = I. 

The excited-state density matrix fi x ^(k) is related to 
thc normal Green's functions via thc relations 



^(k) = limVe-"V 1 e (k, 8Wn ) 

77 — s-0 £ — * 



(8) 



The corresponding lower- and upper-polariton densities 
nip(k) and fi up (\f) are obtained from ([8]) in terms of the 
generalized Hopficld coefficients. Thus, for a fixed to- 
tal polariton density n p , the condensed lower-polariton 
population is given by 



|$i P | 2 = rip - nip 



(9) 



From this quantity, $ x and <I> C are finally obtained via 
the Hopfield coefficients for the condensate Xq and Co , as 
derived from Eq. (J3J) . Hence, a fully self-consistent solu- 
tion can be obtained by solving iteratively Eqs. (J3|), j5]), 
©, and until convergence of the chemical potential 
p, and of the density matrix n X £ (k) is reached. From this 
solution, for a given polariton density n p and tempera- 
ture T, we obtain the spectrum of collective excitations 
j^ip(up) an( j ^ e polariton population ni p r up \(k). 

For the numerical calculations, we use parameters de- 
scribing the experiment in Ref. [2j, i.e. hfl R = 13 meV 
and detuning 6 = e% — eg = 5 meV. We study the results 
as a function of the system area A, ranging from 100 /im 2 
to 1 cm 2 . Thc momentum-dependent interaction poten- 
tials v x and v s 17J, are evaluated for a CdTe quantum 
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well. In Fig. [T]we show the energy-momentum dispersion 
of the collective excitations, ±£^. p arLC1 ±E%. P , as obtained 
for two different values of the total polariton density n p 
above the condensation threshold, at T = 20K. Close to 
zero momentum (see inset), the dispersion of the lower 
polariton branch becomes linear, giving rise to phonon- 
like Bogolubov modes 0, [2(| , as in the standard single- 
field theory [22||. Due to the non-linear terms H x and 
H s , the polariton splitting decreases for increasing n p . 

The interplay between exciton saturation and interac- 
tions in determining the energy shifts has no counterpart 
in the BEC of a single Bosc gas. Here, the two effects 
produce independent energy shifts of the two polaritons. 
We plot in Fig. [2] (a) the energy shifts of the two polari- 
ton modes at k = 0, as a function of the density. Exciton 
saturation and interactions result in a blue-shift of the 
lower polariton and a red-shift of the upper polariton. 
The shifts are linear in n p with the slope changing by a 
factor of two across the threshold (see inset), because the 
contribution of the thermal populations n xx , n xc in Eqs. 
§5§ is twice that of the condensed ones n xx , n° xc . This 
trend and the magnitude of the shifts reproduce fairly 
well the experimental data 0. To explain the origin of 
the opposite shifts of the two polaritons, we plot in Fig. 
[2] (b) the exciton energy Eq = eg + Sff , and the exciton- 
photon coupling Eff, as a function of n p . The two quan- 
tities contribute comparably to the deviations from the 
ideal Bose-gas picture. We predict a very small reduc- 
tion of the polariton splitting up to the largest polariton 
density estimated from the experiments, thus support- 
ing the idea that polaritons - as hybrid exciton-photon 
quasiparticles - are stable well above the BEC threshold. 

We now turn to study the thermodynamical proper- 
ties of polariton BEC. Fig. [3ja) shows the polariton 



fluctuations depleting the condensate in favor of excita- 
tions. 



population ni p (E) 



7 (E) for two values of the to- 



tal density n p , below and above threshold respectively. 
Below threshold, polaritons follow a Maxwell-Boltzmann 
distribution. Above threshold, the distribution becomes 
highly degenerate, with a macroscopic occupation of 
the lowest-energy state, and a saturation of the pop- 
ulation at high-energy. Fig. [3] (b) displays the sim- 
ulated one-body spatial correlation function g^(r) = 
(^( r )v; c (0))/(n cc (r)n cc (0)) 1 / 2 of the photon field. This 
quantity is directly related, via the photon fraction, to 
the actual polariton correlation function, and models 
the outcome of an optical experiment. It shows the oc- 
currence of ODLRO above threshold which is the main 
feature of BEC of an interacting Bose gas 22j. Below 



threshold, the correlation extends only over the thermal 
wavelength At — 1 fJXa.- In experiments the corre- 
lation pattern is shaped by the interface disorder. How- 
ever, the measured long-range correlation is always below 
40%, as compared to 80% of our prediction. By means 
of a kinetic model, we have recently suggested [14J that 
this discrepancy is the main result of deviations from 
thermodynamical equilibrium, with enhanced quantum 




FIG. 1: The dispersion of the normal modes of the system for 
polariton density n p = I5um~ 2 (solid) and n p = 400/xm -2 
(dashed). The non-interacting photon (dash-dotted) and ex- 
citon (dotted) modes are also shown. Inset: detail of the 
low-energy region, showing the onset of the linear Bogolubov 
dispersion (the blue straight line is a guide to the eye). 
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FIG. 2: (a) Lower (solid) and upper (dashed) polariton ener- 
gies at k — vs polariton density n p . Inset: Double log- 
arythmic plot of the lower polariton energy (dotted lines: 
linear slopes below and above threshold), (b) Bare exciton 
energy Eg (solid) and effective exciton-photon coupling Eji 
(dashed). All quantities were computed for T — 20 K. 

In Fig. HJa), we report the computed density- 
temperature phase diagram. In the plot, some values 
of the exciton fraction \Xq\ 2 in the polariton condensate 
are indicated along the phase boundary. It decreases for 
increasing density, due to interactions, but stays finite, 
confirming the stability of polaritons up to high density. 
Fig. 0Jb) shows a detail of the low-T region of the phase 
diagram, computed for different system areas A. We 
also display the physe boundary of the normal-supcrfluid 
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FIG. 3: (a) Polariton population vs. energy for two values 
of n p , computed at T = 20K. (b) Corresponding one-body 
spatial correlation function of the photon field. 
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FIG. 4: (a). Phase diagram of polariton BEC, computed for 
the parameters of Ref. 0], and A = 100 fim 2 . The exciton 
fraction in the condensate |^o| 2 , along the phase boundary, 
is indicated in boxes, (b) Detail of the low-T region. HFP 
denotes the result of the present theory. MF is the mean-field 
result. QC denotes the quasi-condensate transition, corre- 
sponding to the onset of a superfiuid density. 



(quasicondensate) transition, as obtained from an exten- 
sion of the Landau formula [12] . For T < 20 K and 
A < 1000 /mi 2 , the BEC phase boundary lies well below 
the quasicondensate one. Our BEC picture is thus well 
suited for the description of recent samples @, Q , char- 



acterized by polariton localization. In the case of a more 
extended, spatially homogeneous system, a description 
in terms of the Bcrczinski-Kostcrlitz-Thouless transition 
would be necessary. Fig. H£b) also shows the logarithmic 
increase of the critical density as a function of A, due 
to the increase of thermal fluctuations. Quantitatively, 
the variation is very small, in particular for T > 20 K 
■ This dependence and the quasicondensate behaviour 
are hence only expected in samples with improved inter- 
face quality and at lower temperature. The thin lines in 
Fig. 2] (b) are the result of a mean-field approximation, 
obtained by neglecting the excited states population, i.e. 
n ?p + n up — n p- This approximation overestimates the 
group velocity at k = 0, thus strongly underestimating 
the critical n p at low T. This points out to the impor- 
tance of the HFP approach that we have adopted. 

In conclusion, we have generalized the HFP theory to 
the case of two coupled Bose fields at thermal equilib- 
rium. The theory allows modeling the BEC of polaritons 
in semiconductor microcavities in very close analogy with 
the BEC of a weakly interacting gas. The predicted crit- 
ical density is in good agreement with a recent measure- 
ment Q. Our analysis thus supports the interpretation 
of the experimental findings in terms of a transition to a 
quantum degenerate Bose fluid. Open questions remain, 
basically related to the role of disorder and localization. 
If the sample quality can be improved, polaritons will be- 
come an invaluable tool for studying the effects of dimen- 
sionality and fluctuations in interacting Bose systems. 
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